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Abstract 

Ah initio calculations of the structure and electronic density of states (DOS) of the perfect core- 
shell Ag27Cuj nanoalloy attest to its D^h symmetry and confirm that it has only 6 non-equivalent 
(2 Cu and 4 Ag) atoms. Analysis of bond-length, average formation energy, heat of formation of 
Ag27Cuj and LI2 Ag — Cu alloys provide an explanation for the relative stability of the former 
with respect to the other nanoalloys in the same family. The HOMO-LUMO gap is found to be 
0.77 eV, in agreement with previous results. Analysis of the DOS of Ag27Cu'j, LI2 Ag — Cu alloys 
and related systems provides insight into the effects of low coordination, contraction/expansion 
and the presence of foreign atoms on the DOS of Cu and Ag. While some characteristics of the 
DOS are reminiscent of those of the phonon-stable LI2 Ag — Cu alloys, the Cu and Ag states 
hybridize significantly in Ag27Cui, compensating the d-band narrowing that each atom undergoes 
and hindering the dip in the DOS found in the bulk alloys. Charge density plots of Ag27Cu-j 
provide further insights into the relative strengths of the various interatomic bonds. Our results 
for the electronic and geometric structure of this nanoalloy can be explained in terms of length 
and strength hierarchies of the bonds, which may have implications also for the stability of alloy 
in any phase or size. 
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I. INTRODUCTION 



Small bimetallic nanoclusters often have physical and chemical properties that are distinct 
from that of their pure bulk counterparts and suggestive of novel applications. [H 121 E] Not 
surprisingly, materials assembled from finite-sized bimetallic clusters have been investigated 
intensively not only for their catalytic and optical properties, [H El El d |Hl El HDl El 
[T2I [T3l [H] but also for their ability to assemble into duster crystals [151 [IS] and their 
possible applications in single-electron tunnelling devices. |lj Along with its high symmetry 
and relatively high melting temperature, one of the criteria for a cluster to be used as a 
potential building block for cluster-assembled materials is its chemical stability relative to 
other reagents and to other clusters of the same material. Also, major difficulties arise 
from the fact that clusters may tend to coalesce when assembled. This can be prevented 
in one of two ways - either by isolating the clusters in matrices or by coating them with 
surfactants. [TB] An alternative route is to find nanoclusters that are naturally stable, i.e., 
nanoclusters whose intra-cluster interaction is stronger than the inter-cluster interaction 
allowing the clusters to keep their individual identity intact upon assembling. Even so, 
cluster-assembled materials could still be metastable against dissociation into their bulk 
phases. 

Darby et al. [3j, using many-body Gupta potentials, studied the structure and stability 
(as refiected by the total energy) of a wide variety of Cu^Auy nanoclusters with up to 56 
atoms and x/y = 1,3; corresponding to the well-known x : y ratios that result in stable 
ordered bulk phases at low-temperatures. jT^ They found that the geometry of the cluster 
is infiuenced by the tendency to maximize the number of Cu — Au and Au — Au bonds. 
Rossi et al. [Ij, on the other hand, proposed a new family of 34-atom bimetallic alloys using 
the genetic global optimization technique (GGO). fT] These nanoalloys are characterized 
by a perfect core-shell structure in which the smaller atoms {Cu or Ni) compose the core 
whereas, the relatively larger, Ag atoms lie on the surface. They find the nanoalloys to 
be energetically and thermodynamically more stable than pure clusters Ni^/^, Cu^^, and 
Ag34, [1] and attribute the relative stability of the nanoalloy structures to the supplanting 
of the inner Ag atoms by smaller atoms {Cu or Ni) thereby reducing the internal strain in 
^4(734, or the replacement of outer Cu atoms by larger atoms {Ag) to reduce the external 
strain in CU34. As we shall see later, the most stable of the 34-atom Ag — Cu nanoalloy 
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family proposed by Rossi et al, [T] Ag27CuT, provides a hint that it is not the x : y ratio 
that guarantees the stabihty of either bulk or nano-alloys, rather it is the maximization of 
the number of optimized Cu — Cu and Cu — Ag bonds. 

Among a set of possible core-shell nanoclusters modelled by many-body interatomic po- 
tentials, Rossi et al. chose the compositions corresponding to their most stable structures 
for some selected sizes and locally optimized the structures using density functional the- 
ory (DFT) to confirm the trends given by the GGO and to single out the clusters with 
high electronic stability, namely, with the largest width of HOMO-LUMO (Highest Occu- 
pied Molecular Orbital - Lowest Unoccupied Molecular Orbital) gap. They determined the 
thermodynamic stability of the chosen structures through calculations of the melting tem- 
peratures from molecular dynamics simulations and of temperature-dependent probabilities 
of the global minima by harmonic thermodynamics. [18] Among the 34-atom family of nan- 
oclusters, they found Ag27Cu7 and Ag2TNi-j to have the least excess energy with respect to 
bulk atoms (the lowest heat of formation), strong electronic stability (large HOMO-LUMO 
gap), and relatively high melting temperatures. [I] For these compositions, the GGO method 
finds the lowest energy structure to have D^h symmetry in which the 7 coie-Cu atoms form a 
decahedron while the 27 shell- atoms are placed in an anti-Mackay overlayer. In general, 
the formation of nanoclusters is driven by the tendency to minimize dangling bonds. [I] 
Even though such a tendency would lead to spherical structures, [1] the presence of ci-states 
increase the tendency to create directional bonding {s-d hybridization), which is manifested 
by the preference for well-ordered atomic structures, such as the decahedral forms. [TU] 

Experiments and heat of formation calculations have shown that Ag — Cu alloys generally 
tend to segregate. [201 EB [22] In a sense, the core-shell structure of Ag27Cu7 nanoalloy is 
itself segregated. In order to understand how its particular geometry implicitly stabilizes it, a 
detailed examination of the relative strengths and lengths of the Ag — Cu^ Cu — Cu, and Ag — 
Ag bonds is needed. For bulk Cu — Au alloys the presence of a dip in the electronic DOS near 
the Fermi level was also considered to be a signature of alloy stability. [20] Interestingly, such 
a dip is not found in bulk Au — Ag alloys. [TTl |20] Moreover, in the case of crystalline solids, 
the structural stability is linked to the absence of phonon instabilities. [23] The purpose 
of this paper is to carry out a full investigation of the relationship between the geometric 
and electronic structure of Ag2'jCu-j nanoalloy and related bulk systems to get insights into 
the various factors that may impact their stability. That is, through examination of the 
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formation energy, the density of states near the Fermi level, the HOMO-LUMO gap, and 
the charge density distribution, we develop criteria which may lead to chemical and electronic 
stability of Ag2-jCu-j. To obtain additional insights into the structure-stability relationship, 
we have also carried out calculations of the bond-length, the electronic structure, and the 
phonon dispersion of Ag^Cu and Cu^Ag bulk alloys in their LI2 phase. |2l] 

The rest of the paper is organized as follows: Section |TT] contains the computational 



details, while Section |III| is a summary of our results and is divided in subsections |III A 



and IIIB[ Subsection III A analyzes the geometry and bond coordination of the atoms in 
the Ag27Cu-j nanoalloy and those in the bulk Ag — Cu LI2 alloys, for insights into proposed 



stability criteria. In Subsection |IIIB[ we examine our calculated DOS of Ag27Cuj and of the 
bulk systems {CuzAg and Ag^Cu ), and the charge density distribution of Ag27Cui. Finally, 



in Section |IV] we summarize our conclusions, and discuss how the relation between structure 
and stability in Ag2TCuj nanoalloy can be understood in terms of a specific hierarchy in 
bond strength and the capability of the structure to provide the bond lengths for which that 
hierarchy is satisfied. 



II. COMPUTATIONAL DETAILS 



Periodic super-cell calculations are performed in the framework of density functional 
theory. [25] Our calculations are based on the pseudopotential approach and the plane wave 
method (Quantum ESPRESSO: opEn-Source Package for Research in Electronic Structure, 
Simulation, and Optimization). [26] Ultra-soft pseudo-potentials [27] used here are generated 
consistently with GGA schemes. For the GGA functional the expression introduced by 
Perdew, Burke, and Ernzerhof (PBE) has been introduced. [2H] Integrations up to the Fermi 
surface are performed by using a broadening technique [2S] with smearing parameter of 0.2 
eV (0.147 Ry). Below we provide some specifics of the calculations as relevant to a particular 
system. 



A. Calculation of bulk systems 

We have performed extensive convergence tests for lattice constants, bulk moduli and 
total energies of bulk Cm, Ag, Ag^Cu , and Cu^Ag . To obtain the minimum energy config- 
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uration with zero stress based on total energy differences of 1 mRy, as well as, convergence 
up to the third and second digit in the lattice parameter (in a.u.), and the bulk modulus 
(in Mbar), respectively, the calculations demand a k-point sampling of 145 Monkhorst-Pack 
special points [30] (corresponding to a 16 x 16 x 16 Monkhorst-Pack grid) for the integrations 
over the Brillouin zone (BZ). Furthermore, the plane wave kinetic energy cut-off, Ecut, and 
the energy at which the charge density Fourier expansion is truncated, Ep, had to be set 
equal to 680 eV (50Ry) and 8160 eV (600 Ry), respectively. These convergence criteria 
surpass by far most of those reported in the hterature, but are necessary to obtain reliable 
results in the present case, as we will see. 

The lattice dynamics of Ag^Cu and Cu^Ag bulk alloy at arbitrary wave-vectors is ob- 
tained by the Perturbational DFT (DFPT) which is based on the linear response the- 
ory. [311 [321 [33] To determine the force constants we use a 2 x 2 x 2 q-point mesh in 
the BZ of the LI2 structure. [32] Phonon dispersion curves are obtained by the standard 
Fourier interpolation method. [32] 

B. Calculation of Ag27Cuj nanoalloy and isolated atoms 

Since in the unrelaxed configuration of Ag27Cu-j the separation between the most distant 
atoms is about 8.7 A, we locate the nanoalloy inside a cubic super-cell with side length of 24 
A. In this manner we ensure that as a result of periodic boundary conditions, the atoms at 
the edges of neighboring clusters are at least 15 A apart, thereby isolating the clusters from 
each other. The same cubic box is used to model and calculate the total energy of isolated 
Cu and Ag atoms, using a spin-polarized calculation. 

In this work the Broyden-Fletcher-Goldfarb-Shanno (BFGS) algorithm ^] is used to 
minimize the nanoalloy total energy as a function of atomic positions. [25] At equilibrium, 
forces on the nanocluster atoms are required to be below 6.5x10^'* eV/A (2.6xl0~^ Ry/au). 
Given the large separation between nanoclusters (isolated atoms) in neighboring supercells, 
integrations over the BZ using only 1 k-point lead to reliable convergence in the calculated 
values of the total energy (see Ref. For the nanoalloy, E^ut and Ep are 680 eV and 
8160 eV, respectively, as mentioned above. Since these parameters surpass the demands for 
convergence in other systems involving copper or silver, and using the same DFT code, [36l 
l37] we expect them to work well for the nanoalloys of interest here. While these demands for 
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convergence make the calculations very cpu intensive, it is worth mentioning that we find 
the total energy and even the geometry of the nanoalloy to be severely affected if we were to 
use the default values of Ecut = 340 eV and Ep = 1360 eV in the code. For example, with 
the latter choice, the fully relaxed D^h structure (resulting from Ecut = 680 eV and Ep = 
8160 eV) is no longer stable and relaxes towards a structure in which the atomic positions 
break the D^h symmetry by displacements of up to 1 A from their original positions and 
forces cannot be lowered below 2.5 x 10^^ eV/A. 



III. RESULTS AND DISCUSSION 



Since the nanoalloy presents several opportunities for comparison of its properties and 
development of criteria for its stability, we find it beneficial to divide this section into two 



subsections, each of which consists of several parts. In subsection III A, we concentrate on 
issues related to the geometry, the distribution of bond lengths, the atomic coordination, and 
the formation energy of the nanoalloy. We first introduce a notation in subsection III A 1 that 



classifies the atoms in Ag27Cuj according to their location within the nanoalloy. In III A 2 



we inspect how the local coordination of the atoms in the nanoalloy relates to the bond 
lengths. Since there is hardly any experimental data on Ag2-jCu-j and since calculations 
of the phonon density of states of the nanoalloy from first principles is still a challenge, we 



have included in subsection |III A3| our results for the structure and dynamics of bulk alloys, 
Ag^Cu and Cu^Ag, to gain insights and draw stability criteria for the nanoalloy of interest 
here. The average formation energy of Ag27Cuj nanoalloy and stability considerations are 



analyzed in subsection III A 4 



In subsection IIIB we focus on the electronic DOS and the local charge density of the 



nanoalloy. Subsection III B 1 contains the electronic DOS of the Ag2-jCu-j nanoalloy and 
includes for comparison also those of bulk alloys, Ag^Cu and Cu^Ag. The local charge 



density distribution in the nanoalloy is summarized in subsection III B 2 
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A. Geometry, bond coordination, and stability considerations 



1. Geometric Structure of Ag27Cui nanoalloy 

The initial configuration we adopted for Ag27Cuj [38] nanoalloy nicely relaxes towards 
the D^h symmetry after energy minimization, as seen from the plots in Fig. [Tj Accordingly, 
there are only 6 types of non-equivalent atoms: 4 types of Ag and 2 types of Cu atoms. 
This leads to a natural and useful classification of the atoms that refers to their distance 
from the mirror plane: Cu layer (CmO), Ag layer {AgQ), Cu layer ±1 (CmI), Ag layer 
±2 (A(yf2), Ag layer ±3 {Ag?))^ Ag layer ±4 [AgA)^ as shown in Fig. [l] Layer 0, which lies on 
the mirror plane, consists of two pentagonal structures; the smaller one is made of copper 
atoms (CmO) and fits in the larger one that is composed of silver atoms {AgQ) (Fig. [l]^a) and 
(b)). The single-atom layers^ layers ±1 and ±4, sit on the 5-fold rotation axis (see Fig. [l](c) 
and (f)). The other two pentagonal layers, ±2 and ±3, are centered at the 5- fold symmetry 
axis (Fig. [l](d) and (e)). The radii of Ag2 and CuO pentagons are parallel to each other, but 
rotated 36° with respect to the AgQ and Ag2i pentagons, as shown in Fig. |2} Since the layers 
are symmetric with respect to a mirror plane, the cluster can be characterized by only five 
of them, say, the central layer (layer 0) and those above this (layers 1 to 4). Ultimately, the 
symmetry of Ag27Cu7 allows us to fully describe its geometric structure by 8 parameters: 
the interlayer distances (rfoi = 1.341, di2 = 0.052, ^23 = 0.992, ^34 = 1.512 A) and the 
pentagons, side length {acuo = 2.584, a^go = 4.897, aAg2 = 5.115, aAg3= 2.948 A, as shown 
in Fig. [2]). The dislocations in the D^h structure, as present in our initial configuration, [38] 
relaxed into the perfect D^h structure under the stringent criterion that the apothem of each 
pentagon can be well defined up to 0.0003 A, while interlayer distances are well defined up 
to 0.0001 A. The interatomic bond lengths for the 34 atoms in the cluster, in the relaxed 
geometry, are summarized in Table |T} and discussed in detail below. 

2. Neighbor distances in Ag27Cu-j: a comparison with bulk Ag and bulk Cu values 

Despite the perfect D^h symmetry of the cluster. Table |T] shows an intricate hierarchy of 
bond lengths in the optimized structure. For example, Ag3 atoms have 9 neighbors spread 
in a range from 2.726 to 2.948 A, all within a separation equal to the bond length of Ag 
bulk. Figs. [3] (a) and (b) show a comparison between the bond coordination of atoms in the 
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Ag2'jCui nanoalloy and that of atoms in bulk Ag and bulk Cu. Notice from Table |T] that 
the local coordination of Cu atoms in Ag27Cui seems, in fact, not dramatically different 
from that in bulk Cu regarding the number of first and second NN. The CuO atoms, for 
example, have 12 neighbors within a distance of 2.548 and 2.786 A, and another twelve 
between 4.181 and 4.482 A. Silver atoms, on the other hand, find themselves in unusual 
atomic environments: at distances 3 A (bulk Ag nearest neighbors distance) Ag2 and 
Ag4: atoms have barely acquired 6 neighbors, while Ag3 and AgO atoms get 9 and 8 NN, 
respectively. Between 4.2-4.5 A, at which bulk Ag atoms already have 18 NN, AgO, Ag3 
and Ag4 atoms have acquired at most 11 neighbors, while Ag2 atoms have only 7 NN, 
finding themselves as the most under-coordinated atoms of the cluster. Most importantly, 
we notice that the NN bond lengths between Cu atoms contract by as much as 2% with 
respect to the value in bulk (2.599 A). Such contraction may be expected for Cu atoms in 
low coordination environments, such as on the surface layer of Cm(IOO) (inward contraction 
of around 3 % |23]), and not in Ag2jCu'j because of their high coordination (12) as in bulk 
Cu. Interestingly, for the shell Ag atoms, which have much lower local coordination than 
the core Cu atoms, the Ag — Ag NN bond lengths are at most ~ 2.6 % shorter than those in 
bulk Ag (2.943 A). If coordination alone were a measure, one might have expected the most 
under-coordinated Ag atoms, Ag2 and AgA, to undergo a larger contraction, as found on Ag 
surfaces. [23] Instead, in Ag2jCuj, Ag2 and Ag4 form short bonds of about 2.59 A with their 
neighboring Cu atoms, as seen from Table |Tj These bond lengths are in fact very close to the 
smallest Cu — Cu bond lengths (2.55 A) and considerably smaller than the smallest Ag — Ag 
bond lengths (2.87 A). Indeed, it follows from Table [l] that in Ag27Cuj the first NN of every 
Ag atom is a Cu atom, pointing to the reality that finite sized structures of these elements 
may not follow straightforwardly the behavior of infinite and/or semi- infinite systems and 
the relationship between bond coordination number and bond stiffening might be subtler 
in nanoalloys. Two conclusions may, nevertheless, be drawn from the above: 1) the fact 
that Cu — Cu bond lengths in Ag27Cu7 contract almost as much as Ag — Ag bond lengths 
suggests that Cu atoms are more sensitive than Ag atoms to either low coordination, or 
local geometry, or chemical environment; 2) the low coordination of the Ag atoms appears 
to be significantly compensated by the formation of short bond lengths with Cu atoms. 
To discriminate between the above mentioned effects of local coordination, geometry, and 
environment, the conclusions about bond lengths in the Ag27Cu7 nanoalloy need to be put 
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on firmer grounds through examinations of details of the electronic structure and the charge 
density distribution, and their implications for nanoalloy stability. We will turn to this in 



Sections III A 4 and IIIB But before we do that, it is interesting to examine the structural 
stability of related bulk alloys, Ag^Cu and Cu^Ag, for which some information already 
exists and hence can serve as reference points. The relationship between the composition of 
Ag2iCuj and Ag^Cu is obvious. The other alloy is chosen to establish whether preponderance 
of Cu and/or its effect on bond lengths is a key for understanding the structural stability of 
these alloys. 

3. Structure, phonons, and heat of formation of Cu^Ag and Ag^Cu bulk alloys 

In considerations of structural stability of bulk alloys, it is essential that the heat of 
formation be negative and that the phonon spectrum be well defined. Ag — Cu alloys, unlike 
Au — Cu and Au — Ag alloys, are known for their tendency to segregate and have a miscibility 
gap beyond the eutectic temperature of the material. [201 EH ES] Earlier calculations [20l 
[2T| |22] have shown that Ag — Cu alloys possess a positive heat of formation regardless of the 
chosen ratio of Ag to Cu. [201 El] In particular, the heat of formation per atom was found 
to be about 70 meV for Ag^Cu (60 meV, in this work) and 80 meV for Cu^Ag (64 meV, in 
this work), pointing to the structural instability of these Ag — Cu bulk alloys. |2I] 

It has, however, been pointed out that even with a positive heat of formation presence 
of a well-defined phonon spectrum may serve as an indicator of alloy stability under special 
formation conditions, for example, using non-equilibrium techniques. [391 SOI SIl 111] Indeed, 
among a variety of immiscible noble-transition metal alloys, some - particularly those that 
mix fee and hep metals - have shown mutual solid solubility. [43] Kong et al. ^39j have 
calculated the phonon spectra of several structures of the equilibrium immiscible Ag^Ruy 
alloys, for x/y = 1/3 and 3, to find that only the LI2 and -DO19 phases of Ru^Ag may be 
stable. In Ref. [39], phonon-stable/unstable structures were associated with the presence of 
relatively high charge density bridging NN atoms of the same/different element. We will 



come to charge density implications later in Section III B 2 For the present discussion, 
we note that, regardless of the structure of the stable phases, the phonon-stable Ag — Ru 
alloys are all Ru-iich and have smaller lattice parameter than the corresponding Ag-iich 
structures, owing to the fact that typical bond lengths in bulk Ru are smaller than those in 
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bulk Ag. Moreover, Kong et al. were able to remove the phonon instabilities, when present, 
by artificially increasing the external pressure; i.e., by simply reducing the equilibrium lattice 
parameter. 

In the above spirit, we now turn to the calculation of the structure and the lattice 
dynamics of Ag^Cu and Cu^Ag bulk alloys using DFT and DFPT methods. We find the 
bond-length for Ag^Cu to be 2.87 A and that of Cu^Ag to be 2.70 A. For reference, note 
that our calculated bond lengths for bulk Cu and Ag are 2.59 and 2.94 A, respectively. Note 
also that the shortest Ag — Cu bond lengths in Ag2'jCu^ are around 2.6 A. The bond-lengths 
in bulk alloys Ag^Cu and Cu^Ag are thus larger than the shortest Ag — Cu bond-length 
in the nanoalloy and that in bulk Cu and may imply lack of overlap of the cZ-orbitals for 
its Ag — Cu and Cu — Cu bonds (as we shall see), pointing to the structural instability of 
these bulk alloys. Our calculated phonon dispersion curves (Fig. |4]) of Ag^Cu and Cu^Ag, 
showing the absence of unstable modes, however, suggest that these alloys may be stable 
and obtainable by non-equilibrium techniques. [SH] 

Phonon dispersion curves are furthermore a measure of the contribution of the vibrational 
entropy to the free energy of a given system. [S] The vibrational entropy integrates the 
vibrational DOS weighted by a factor that falls off as the frequency of phonons increases. 
Thus, vibrational entropy plays a larger role in the minimization of the free energy for 
systems whose (well-defined) phonon dispersion curves display notable contributions and 
shifts of the density of states towards the lower frequency range. [22] For example, bulk 
Cu^Au {LI2) has larger vibrational entropy than either bulk Au or Cu. [45J Our calculated 
phonon dispersion for bulk Ag^Cu (see Figs. |4]^a)) is softer than that of bulk alloy Cu^Ag 
(Figs. |4](b)), bulk Ag, ^46j and bulk Cu, |16] indicating that the vibrational entropy of 
Ag^Cu is larger than that of the bulk Ag and Cu. Regardless, the vibrational entropic 
contribution to the reduction of the free energy is small [22| (typically less than 5 meV at 
room temperature) compared to the positive heat of formation (~60 meV) found for this 
bulk alloy. Instead, the presence of stiffer vibrational DOS of Cu^Ag, compared to that of 
Ag^Cu, reflects stronger bonds in the former which may provide stability if created by the 
techniques mentioned above. [39l HOj HH |12] 
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4- Formation energy of Ag27Cuj 



In Ref. pp, the thermodynamic stabihty of a given nanoalloy is evaluated via consider- 
ations of its melting temperature and its relative energetic stability is established through 
comparison of the heats of formation, A, adapted to binary clusters as follows: 



where N = N1 + N2, E[AgN^CuN2) is the total energy of the nanocluster, and E{Cubuik) and 
E{Aghuik) are the total energies of one Cu and one Ag atom, respectively, in the bulk phase. 
Remarkably, although clusters with increasing binding energies per atom do not necessarily 
have higher melting temperatures, [37] Ag2'jCuj came out with both the highest melting 
point and the least heat of formation. In terms of stability and minimum-energy structures, 
however, it is important to know also the output given by the average formation energy, 
which measures the dissociation or cohesive energy of the nanoalloy, and to analyze the 
meaning, implications, and scope of these two energetic considerations. Thus, to estimate 
the average strength of the bonds, we calculate the so-called average formation energy per 
atom, Eform, which is defined as. 



where E{Cufree) and E{Agfree) are the energies of isolated Cu and Ag atoms, respectively. 
We find E form{^927Cu-j) to be 2.17 eV. Since no experimental data on this particular binary 
nanocluster exists, we turn to the formation energy of related systems for comparison. For 
example, our calculated cohesive energy of bulk Ag, Ag^Cu, Cu^Ag, and Cu are 2.51, 2.66, 
3.06, and 3.34 eV, respectively, implying that Eform{Ag27Cu-j) is smaller than all these bulk 
values. Note that the higher cohesive energy of bulk Ag^Cu as compared to that of bulk 
Ag signals a stronger Ag — Cu bond than the usual Ag — Ag one. In fact, from the results 
presented in Ref. [T], one finds that the average formation energy per atom of the 34-atom 
family decreases monotonically from ~2.6 to ~2.0 eV as the Ag content increases from 
to 34. Considerations of formation energy alone would thus imply that in this family of 
nanoalloys the Cu — Cu and Cu — Ag bonds are stronger than the Ag — Ag bonds and 
that Ag27Cuj is not the most stable structure. It is thus surprising that a related quantity, 
i.e. the heat of formation (defined as in Eq. |2| but substituting E{Cufree) and E{Agfree) 



A^AgNiCuN^) 



E{AgN,CuN,) - NiEjAgbuik) - N2E{Cubuik) 

nI 



(1) 



Eform {^dNi C UN2 ) 



E{Ag^,CuN,) - N,E{Agfree) - N2E{CUfree) 

N 



(2) 
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by the cohesive energy of Cu and Ag, respectively), plotted as function of the Ag/Cu ratio 
in Ref., [T] shows a minimum at this intermediate composition {Ag2'rCuT) - a result that 
stands in contrast to that found in Ag — Cu bulk alloys. j^Dl Ell 122] Nevertheless, we find 
that the heat of formation of the Ag27Cu'j is 9 times larger than that of Ag^Cu, but this 
suggestion of instability might be misleading since the heat of formation not only measures 
the strength of the bonds but also weighs the energetic cost (gain) of breaking(forming) 
single element bulk bonds to form(from breaking) binary bonds. In reality, the formation 
of intricately tailored structures as ^(727(7^7 is not expected to occur simply by melting the 
parent compounds. Thus measures like heat of formation have to be supplemented by others 
such as the dynamical stability of the alloy as displayed by its vibrational modes. The heat 
of formation is perhaps more of an indicator of the life-time of the nanoalloy, say, against 
clustering and the eventual formation of segregated metallic bulk Cu and Ag, if one ignores 
the energy barriers needed to actually break all bonds in the nanoalloy. 

On the experimental side, in addition, we note that the pure clusters Agj (Ref. 47) and 
Ag^g (Ref. [1H|) were found to have dissociation energy of 2.73 eV and 2.88 eV, which are very 
close to the experimentally observed cohesive energy of bulk Ag (2.94 eV). The formation 
of pure cluster structures, such as Ag^ and Ag^g, may thus be seen to be energetically 
more favorable than the nanoalloys. Perhaps, the possible disintegration of Ag27Cuj into 
pure-element clusters may be argued against on the basis of the strength of the Ag — Cu 
bond. To estimate the strength of the bonds in the 34-atom nanoalloys and understand 
what distinguishes Ag2-rCuj in its family of nanoalloys, we turn to Fig.2 of Ref. pQ. Rossi 
et al. show that as the amount of Cu increases up to ~20 % (starting from Ag^^), the 
heat of formation is reduced or kept constant, implying that small amounts of Cu atoms 
immersed among Ag atoms (in the nanoalloys Ag^^^nCun as n decreases from 7 to 1) create 
Cu — Ag and Cu — Cu bonds that are stronger than those in bulk Cu and are able not 
only to counterbalance the cost of the cohesive energy of the newly added Cu atom, but 
also to increasingly stabilize the nanoalloy. Fig.2 of Ref. [1] shows also that if the content 
of Cu increases beyond 7 atoms, the heat of formation increases again, with the implication 
that the strength of the bonds is not able to compensate the bulk Cu cohesive energy for 
an additional atom. In conclusion, Ag2-rCu-j possesses the composition and geometry that 
maximizes the number of Cu — Cu and Cu — Ag bonds using the minimum number of Cu 
atoms. 
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The structural stability considerations presented above are derived entirely from the 
energetics of Ag2-jCu-j. The contributions of vibrational entropy could be important [IS] and 
may lead to a minimum of the free energy (as a function of the Ag / Cu ratio) which is different 
from that in the plot of Rossi et al. They, however, argue against such a possibility. [Ij 
Very recent calculations [19] based on many-body interaction potentials also confirm that 
vibrational entropic contributions are small for this set of 34-atom nanoalloys. Calculations 
of the phonon frequencies of these nanoalloys from first principles, as presented in section 
A. 3 for the LI2 bulk alloys, are still desirable, as such a study could serve as an indicator 
of stable compositions. However, the resulting contributions of the vibrational entropy to 
the free energy are expected to be small and should not change the conclusions drawn 
here. Configurational entropy may in general also play a role in determining the stable 
configurations of these nanoalloys, but it is not expected to be important for Ag27Cu-j 
because of the lack of degenerate isomers. [50] 

B. Electronic structure and charge density distribution 

1. Electronic DOS of Ag27Cuj nanoalloy and LI2 Ag — Cu alloys 

In our calculations the Ag27Cuj nanoalloy has 374 valence electrons occupying 187 dis- 
crete, flat, and double degenerated bands, as shown in Fig. [5j As we can see, the discrete 
(i-bands are densely concentrated between -5.0 and -1.5 eV below the Fermi level. Such an 
ideal behavior of the DOS may not be observed in experiments. [19] One of the reasons is 
the Jahn- Teller deformation which lifts the degeneracy of the levels, smears out the ideal 
DOS, and may also lead to the opening of gaps not predicted by the present approach. [T^] 
The second reason is the actual atomic positions occurring in a cluster due to deviations 
from the ideal structure further modifies the DOS. [I9] For example, Ref. [T9] shows that the 
photo-electron spectrum of Cu^^ (among others smaller and larger) is significantly smoother 
than 187 sharp peaks that one may expect, though it does displays smooth dips. To obtain 
the electronic DOS of bulk Cu^Ag and Ag'^Cu alloys and Ag2-jCu-j nanoalloy from ah initio 
electronic structure calculations, the states are thus broadened using Gaussian functions of 
width 0.14 eV. 

Our resulting DOS of Ag27Cuj is shown in Fig. |6](a)-(c), while that of the bulk alloys 
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is presented in Fig. [6[d) and (e). First of all, the HOMO-LUMO gap (A in Fig. [6[a)) of 
Ag2-rCui is found to be 0.77 eV in the ground state, which is only slightly smaller than that 
reported by Rossi et al, fT] 0.82 eV. As expected, the s-states have negligible contributions 
and the displayed structures have mostly (i-character between -5.3 and -1.5 eV. Fig. [6|b) 
shows that the centroid of the nanoalloy valence band red-shifts ~ 1 eV as 

compared to bulk Ag^ and blue-shifts ~ 0.5 eV with respect to bulk Cu. Fig. |6]^c) shows 
that even though the amount of Cu in the Ag2-jCu-j nanoalloy is ~ 4 times less than that 
of Ag^ it contributes to blue-shift the valence band. As shown in Fig. [6]^c)-(e) and in 
Refs. pn 1^ . the role of Cu is, in general, to enhance the DOS at the top of the valence 
band and to blue-shift the centroid, while the effect of Ag is the opposite. Similar results 
have been reported for Au — Pd nanoclusters, [52] in which the increasing content of Au 
on Pd clusters reduces the density of states at the Fermi level. Fig. |6](b) and (c) show in 
addition that the d-band of the nanoalloy is as broad as that of either pure bulk constituents 
- a point worth of noticing since atoms in a low coordinated environment generally exhibit 
a valence band narrowing. [53j Pure Ag clusters and Cu^Au surfaces, [SH ES] for example, 
have shown this effect. The hybridization of Ag and Cu states in Ag27Cur thus compensates 
the d-band narrowing that each atom undergoes. 

In general, the features in the electronic DOS that discriminate stable alloy phases are not 
yet fully understood. Although it is well accepted [20j that the dip in the DOS at the Fermi 
level is related to the stability of a particular alloy phase, there is no obvious correlation 
between the two since stable and ordered Cu — Au alloys present a dip in the DOS while 
Au — Ag alloys do not. [20j In the former the dip in the DOS changes position, width and 
depth with composition, and structure, and has been found to be related to the electronic 
specific heat. [17] Kokko et al. noticed also that the dip is considerably lessened in the 
layered (tetragonal LIq) CuAuI phase (which reduces the Au — Cu bond) with respect to 
Cu^Au and Au^Cu. Also, based on their electronic specific heat calculations, they infer 
that the dip is even smoother in the disordered phases. In this work, we find that the 
nanoalloy Ag27Cuj displays a less pronounced dip (deepest at ~ - 2.25 eV) in the cZ-band 
region (between ~ -5.3 and ~ -1.5 eV) - where Cu and Ag states hybridize - than those 
displayed by LI2 Cu^Ag and Ag^Cu alloys (deepest at ~ -3.6 eV and ~ -2.5 eV, respectively. 
Compare Fig.|6|^a), (d) and (e)), which are similar to, but broader than, those found in stable 
Cu — Au alloys. [17j Yet, from considerations of the heat of formation, the Ag — Cu alloys are 



14 



marked as being immiscible. Furthermore, we find that the DOS of the Ag atoms (Fig. [6[b) 
and (c)) in Ag2iCui resembles to some extent that of bulk Ag, despite being highly under 
coordinated, while that of the fully coordinated Cu atoms is strikingly different from bulk 
Cu presumably, since half of its neighbors are Ag atoms. 

In order to understand the correlation between the Ag/Cu content ratio, the consequent 
decrease/increase of the bond lengths, and the structure of the electronic DOS, we turn to 
examination of the changes experienced by the DOS of individual Ag and Cu atoms in a 
set of environments (some natural, some artificial): Ag27Cui nanoalloy, a free standing Ag 
and Cu monolayer, bulk Ag — Cu alloys, and bulk Ag and Cu (Figs. [7]j9]). For the bulk 
systems, we consider also the effect on the DOS (Fig. [s] of expanding and contracting the 
lattice constant from the equilibrium value found in our DFT calculations. The DOS of Ag 
and Cu atoms in Figs. [7]j9] allow comparison on a one-to-one basis and not as percentile 
contributions as presented in Fig. |6| 

From Figs. [7|^a) and (b), we note that decrease in the content of Ag (and thus the lattice 
parameter) gives rise to a sharp peak at intermediate energies while the bottom of the d- 
band remains almost unchanged and states from the top of the bulk Ag band (Fig. [7](e)) 
retract to lower energies. A few states, however, appear between -1.5 and -2.5 eV, above the 
range of the DOS of bulk Ag, and hybridize with Cu states. As we shall see, the appearance 
of these higher energy states is related purely to the presence of Cu, whereas the depletion 
of the top of the Ag (i-band results from both the presence of Cu and from the enhanced 
overlap of Ag — Ag orbitals as a result of the decrease in the Ag — Ag bond length. 

To isolate the effect of the bond length, we turn to the DOS in Fig. |8] for Ag and Cu 
atoms in bulk environments with bond lengths different from the equilibrium values. In 
order to maintain a reference point, we have taken the lattice constant of expanded Ag^Cu 
and Cu^Ag to be that of bulk Ag, corresponding to an expansion of 2.4 % in the former 
and 8.9 % in the latter. Similarly, to infer the effect of lattice contraction we have used 
the lattice constant of bulk Cu for Ag^Cu (contraction of 9.8 %) and the bond-length of 
CuO — Cul for CusAg (contraction of 5.5 %). By comparing Figs, and|8](f), we find 
that contraction of Ag — Ag bonds pushes the bottom of the Ag d-hand to lower energies 
and depletes the top of the band, while comparison of Figs. |7(a) with gb) and |7|(b) with 
^d) indicates that the contraction of Ag — Cu bonds also pushes the bottom of the Ag 
(i-band to lower energies and significantly lessens the highest features of the DOS of Ag 
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(i-band in bulk alloys. We also conclude from Figs, [sj^b) andjsj^d), and Fig. [Tj^e), [7]^a) and 
that the presence of Cu is responsible for the appearance of Ag states above the top of 
the bulk Ag band hence improving the overlap with the Cu d-hand. On the other hand, as 
seen from Figs. [7]^e) and[8|e), the expansion of Ag — Ag bonds depletes the bottom of the 
DOS and slightly enhances the DOS at the top. Similarly, the expansion of the Ag — Cu 
in the bulk alloys depletes the bottom of the DOS and enhances the highest peaks of the 
DOS (compare [7|^a) with[8|^a) and[7](b) with|8](c)). Atomic low coordination and expansion 
of the bond length cause the same effects on the DOS of Ag, though these are augmented 
in the former situation. Such was found to be the case of a free standing (111) monolayer 
(Fig. [7](d)), in which the 6-coordination of Ag atoms causes a strong depletion of the bottom 
and enhancement at the top of the d-hand. 

From the above, we conclude that the low coordination of Ag atoms in the Ag27Cu-j 
nanoalloy can indeed account for the depletion of the bottom of the d-hand (Fig. [Tj^c)). 
On the other hand, the reduction of the DOS at the top of the rf-band of AgQ (Fig. [9]^a)), 
Ag2) (Fig. |9](c)), and AgA (Fig. ^d)) atoms, suggests that the presence of Cu at such short 
distances outweighs the effect of their low coordination and the DOS at the top of the 
Ag d-band is not enhanced. The DOS of Ag2 is slightly more complex. As mentioned 



in subsection IIIA2, Ag2 is the most undercoordinated atom of the cluster (See Fig. |3]). 
Consistently, it has the least number of states below -4.5 eV and highest DOS above ~ -2.8 
eV, implying that low coordination effects on the DOS of Ag2 are not completely washed 
out through the short bond with CmO, each one of which is incidentally shared by two Ag2 
atoms. At the same time, as seen in bulk alloys, the hybridization of Cu and Ag states is 
improved by setting off the occupation of states above the top of the d-hand of bulk Ag (see 
Fig.^c) and Figs.ga)-(d)). 

We now turn to the issue of the DOS of Cil atoms in the Ag27Cu7 nanoalloy which differ 
prominently from that of bulk Cu. We consider three different aspects that may influence 
their electronic structure: the conspicuously disparate overall geometry, the presence of Ag, 
and the existence of bond lengths longer than that of bulk Cu. To understand each one 
of these we first note that if Cu — Ag bonds are longer than that of bulk Cu, then the 
bottom of the cZ-band of Cu is strongly depleted and an increasingly sharper peak at the 
top is created, which slightly shifts towards higher energies (compare Figs, with |8|^g) 
and[7](g) withjsj^i)). As already noted for the Ag atoms, the effect of expanding Cu — Cu 
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bonds is similar to that of Cu atoms in a low coordinated environment (compare Fig. |8](k) 
and [rj^i) with [7|^j)), that is, the bottom of the Cu d-hand (states below ~3.5 and ~ 4.0 
eV in Figs, ^i) andjsj^k), respectively) is entirely extinguished. Interestingly, the (i-band of 
the free monolayer is sharply localized at the top edge, resembling that of CuO (Fig. |9](e)) 
and Cul (Fig. [9|^f)) atoms notwithstanding their full coordination in Ag27Cu'j. Contracting 
Ag — Cu bonds in bulk alloys reduces the DOS at the top, which is pushed to lower energies, 
and enhances the bottom of the rf-band, thus improving the Ag — Cu hybridization, as seen 
by comparing Figs. [7](f) withjsj^b), and|8](j). Similar features are found by contracting 
the Cu — Cu bond is bulk Cu (compare Figs. [7]^j) and[8|^l)). More importantly, note in 
Figs. [S^j) that even when Cu — Cu and Cu — Ag bonds in bulk Cu^Ag are shorter than 
that of bulk Cu, the region from ~ 4.0 eV . ~ 5.0 eV is nevertheless strongly quenched, 
suggesting that the presence of Ag also strongly depletes most of the bottom of the (i-band of 
Cu atoms (something similar occurs in Au — Cu, see Ref. [IZj). Observe also that presence of 
Ag introduces states below the bottom of the DOS of bulk Cu (Figs. [7||^f)-(h) and[8|^g)-(j)), 
which hybridize with Ag. 

From the above, we conclude first that the sharp peak (Fig. [7](h)) at the top of the DOS 
of Cu atoms in Ag27CuY nanoalloy, characteristic of low-coordination (Figs. [7](i)), can only 
be accounted for by the relatively long distances (~ 2.7-2.8 A) between Cu atoms and half 
of their nearest neighbors - all Ag atoms - (see Figs. |3](a) and Table I), seemingly leading 
to a weak interaction of the Cu atoms with those far-lying Ag neighbors, as occurs for Cu 
atoms in bulk Ag-^Cu and Cu-^Ag alloys. The DOS of Cu atoms in bulk Ag-^Cu (Figs. [Tj^f)) 
thus indicates that the strength of the Ag — Cu bond is considerably weak for Cu atoms; 
in fact, expanding the lattice parameter (Figs. [Sj^g)) changes insignificantly the DOS of Cu. 
Second, the effect of the Ag environment on the DOS of Cu atoms in Ag27Cui is connected 
with the depletion of states between 2.5 and 3.5 below the Fermi level, i.e., the dip around 
3.0 eV in Fig[7]^h), ^e) and (f). We note in addition that the DOS of Ag27Cuj below the 
dip (~ 3.5 eV) is remarkably high, as compared to that of bulk alloys (see Figs. [Tf^f) and 
(g)) - despite the low Cu content - and generates a much stronger hybridization between Cu 
and Ag states, contrasting that observed in bulk alloys (see Figs. |6](c)-(e)); interestingly, the 
DOS of states of Cu atoms at 3.5 eV below the Fermi level in the compressed lattice of bulk 
Ag^Cu (Figs, jsj^b)) is almost as high as that of Cu atoms in Ag2'jCu'j and also results in 
a stronger hybridization between Cu and Ag states compared to that found in equilibrium 
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bulk Ag^Cu bulk. In the third place, we conclude that the optimum Cu — Ag bonding in 
Ag2-rCui comes about in terms of the electronic DOS through the shortening of Ag — Cu 
bonds that allows the hybridization of the rf-states of Cu and Ag atoms, suppressing by this 
means the dip featuring in the DOS of Cu^Ag and Ag^Cu alloys. 

Finally, we remark that although the DOS of Ag atoms is not changed as drastically as 
that of Cu atoms in bulk Ag — Cu alloys - suggesting that Cu is more sensitive than Ag to 
the chemical environment -, the vulnerability of Cu to the presence of Ag is intermixed with 
long bond-length effects. That is, the presence of Cu depletes the top of the Ag d-hand as 
much as Ag depletes the bottom of Cu d-band. However, the short Ag — Ag and Ag — Cu 
bonds (with respect to bulk Ag) as present in both bulk alloys broaden the rf-band of Ag 
atoms, compensating in this manner the effect of the Cu atoms. On the other hand, the 
Cu — Cu and Cu — Ag bond lengths induced in both bulk alloys never become shorter than 
that of bulk Cu, rather it is quite the opposite. As a result, narrowing of the d-band of Cu 
atoms is triggered, aggravating the depletion caused by the presence of Ag atoms and, hence, 
exaggerating the actual chemical susceptibility of Cu to the Ag environment. In summary, 
the electronic DOS is found unambiguously related to the bond lengths held in a particular 
geometry. 

If we were to extrapolate the above results to related systems, we would speculate that 
the dip in the DOS identifies less stable phases of noble metal alloys. For example, based 
on the patterns found in Ag27Cui and the LI2 bulk Ag — Cu alloys, we predict that the 
stability Au — Ag alloys and the absence of a dip in their DOS [2U] is attributed to the fact 
that the lattice parameter of bulk Au and Ag are nearly identical and that the rf-band of 
bulk Ag lies within that of bulk Au, assuring significant (i-band hybridization. Likewise, 
the hybridization between Cu and Au states is strong in bulk Au — Cu alloys because the 
(i-band of Cu also lies within that of bulk Au, albeit in the region near its Fermi level, while 
that of Ag lies deeper. The overlap between the bands of Ag and Cu, on the other hand, 
is relatively small resulting in weaker Cu — Ag hybridization in bulk Cu — Ag alloys. The 
dip in the DOS of Au — Cu is probably not a sign of stability, rather it may be a sign that 
structures and/or compositions allowing shorter Au — Cu would be more stable (amorphous 
phases, perhaps). The heat of formation of bulk Cu-^Au, for example, is negative because 
the strength of the Au — Cu bond is larger than that of bulk Cu and Au, even at the 
distances dictated by the LI2 phases, [21] but which are not necessarily the distances at 
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which the strength of all three Au — Cu, Au — Au, and Cu — Cu bonds is optimized with 
the corresponding hierarchal importance, as exemplified by LI2 bulk Ag — Cu alloys whose 
long Ag — Cu bonds contrast the significantly short ones found in the Ag27Cu'r nanoalloy. 



2. Charge density distribution of Ag27Cu7 nanoalloy 



The first aspect that comes to mind in Fig. 10 15 of the plotted charge densities of the 



nanoalloy Ag27Cu-j is that Ag atoms barely supply charge to the surface of the nanoalloy. The 
question is whether the surface charge depletion coincides with charge redistribution from Ag 
atoms to Cn atoms since two neighboring metals with significantly different work function 
can give rise to electron transfer from one metal to the other, as reported in calculations 
of Pd clusters on Au{lll) by Sanchez et al, [52j in correspondence with a work function 
difference of ~ 0.3 eV. Here, the work function of Cu is larger than that of Ag by ~ 0.2 



eV. Indeed, Figs. 10 and 11 display higher charge density around Cu atoms than around Ag 
atoms. 

The plot of the charge density in a plane passing through 10 Ag and Cu atoms in Fig. 10 
illustrates that the Cu0-Ag2 and Cul-AgA bonds are linked by the highest bonding charge 
density, corresponding indeed to the shortest Ag — Cu bond lengths (~ 2.58 A), and implying 
that they are stronger than the Cul-Cul ( 2.58 A), CuO-CuO (2.58 A) bonds and even the 
CuO-Cul bonds (2.55 A), in that order (see Figs 



10 



and 



11). The next strength of bonding 



charge density occurs for CuO-AgO bonds, followed by Cul-Ag3 bonds, as shown in Figs. 10 



and 11_ In this case, the bond lengths are 2.72 and 2.73 A, respectively. Notice that, in the 
second set of bonds, the charge density is considerably lower than in Cu0-Ag2 and Cul-AgA 
bonds, which supports the assumption that Ag — Cu interactions die out very rapidly (see 



Section III B 1 ). The third place in bonding charge density corresponds to Cu0-Ag3 (Fig. 12 ) 



bonds and the shortest Ag — Ag bonds: Ag0-Ag3 and Ag2-Ag2 bonds (Fig. 10) whose bond 
lengths are 2.79, 2.87 and 2.89 A, respectively (Table The first might infiuence very little 
CuO atoms since the bond is quite large. However, this bonding charge density appears as 
large as that of Ag0-Ag3, which is 0.1 A further apart. The next in bonding charge density is 



the Ag2-Ag3 bond (Fig.[l3|, whose length is 2.90 A, followed by the Ag3-Ag4: bond (Fig.jToj), 
whose length is 2.93 A. The latter bond length is close to that of the Ag0-Ag3 bond (and 
almost identical to the bond length of bulk Ag); however, the charge density bridging these 
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atoms is slightly weaker. Note that Ag3-Ag3 bonds (Fig. 14) are only slightly longer (2.95 
A) but the bonding charge density around it is less than that around the AgS-AgA bonds. 



The charge density around the Ag0-Ag2 bond (Fig. 15) is substantially lower than that 
around Ag3-Ag3, consistent with a larger bond length, 3.0 A. The next larger bond lengths 
are more than ~ 4 A, which are expected to provide much weaker direct interactions that 
will not be discussed here. 

In reference to the importance of bond strength and length hierarchies in alloys, men- 



tioned in subsection III B 1 , we turn to the discussion in III A 3 of Ref. [32] in which Ag — Ru 
phonon-stable structures were associated with a high charge density bridging atoms of the 
same element, whereas phono n-unstable structures were associated with a high charge den- 
sity bridging atoms of different elements. Instead, following the arguments above, we pro- 
pose that their charge density plots indicate that Ag-nch structures do not allow for strong 
Ru—Ru bonds and so, although strong Ru—Ag and Ag — Ag bonds are present, the structure 
is unstable. Additional support to our assertion comes from the fact that phonon instabilities 
are not present in some Ru-iich structures, [22] while those in Ag-iich structures disappear 
by reducing the equilibrium lattice parameter. [39] Besides, Ag-iich. structures with reduced 
lattice parameter show clearly that the charge density bridging Ru — Ru bonds is enhanced 
(see Fig. [3]^e) of Ref. [39]). The above suggests that in Ag — Ru alloys the strength hierarchy 
is in the following order: {Ru — Ru, Ag — Ru, Ag — Ag}. In fact, relatively large differences 
between the strength of Ru — Ru and Ag — Ru bonds and a pronounced deep dip in the 
DOS can be expected to explain all their results simultaneously. 



IV. SUMMARY 



We have presented DFT calculations of the Ag2'jCu-j nanoalloy to understand its structure 
and relative stability via considerations of its energetics, electronic DOS, and charge density 
distribution. The local coordination of Cu atoms is similar to that in bulk Cu regarding 
the number of first and second NN, whereas Ag atoms find themselves in a low-coordinated 
environment but, in exchange, form Ag — Cu bonds which are short as the shortest Cu — Cu 
bonds. On the other hand, the electronic structure of the Cu atoms in Ag2jCuj deviates 
much more from that of atoms in bulk Cu, as compared to the corresponding case of Ag 
atoms in this nanoalloy. 
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Related bulk alloys, Ag^Cu and Cu^Ag, are found to have positive heat of formation 
and form larger bonds than the shortest ones found in the Ag27Cuj. However, we find that 
the resulting interatomic bonds in these bulk alloys are sufficiently strong such that their 
cohesive energy is larger than that of bulk Ag and their phonon dispersion curves do not 
display instabilities. 

From our analysis of the geometric and electronic structure, we conclude that the relative 
stability of Ag^-jCuT, among its nanoalloy family, is the result of the maximization of the 
number of Cu — Cu and Cu — Ag bonds, using the minimum number of Cu atoms. The 
core-shell Ag — Cu nanoalloys do not behave differently from the corresponding bulk alloys 
regarding segregation tendencies and migration of Ag to the surface, as pointed earlier, [M] 
since Ag27Cu-j is segregated by construction. Furthermore, the segregated structure is the 
attribute that leads to its relative stability, provided the core-shell structure allows formation 
of strong Cu — Cu and Cu — Ag bonds without contracting the typical Ag — Ag bond or 
(most importantly) stretching the typical Cu — Cu bond. 

The HOMO-LUMO gap of is found to be 0.77 eV. The DOS of Ag2^Cu^ shows 

features similar to those of Ag — Cu bulk alloys. We find that the novel features of the 
DOS of fully coordinated Cu atoms in Ag27Cu-j are caused by the relatively long distance 
separating Cu atoms from half of its first NN. Short Ag — Cu bond-lengths, on the other 
hand, improve the hybridization of Cu and Ag states; explaining why the hybridization 
of Cu and Ag states in the Ag27Cuj nanoalloy is stronger than in Ag — Cu bulk alloys. 
The observed differences in electronic DOS between Ag27Cui and LI2 alloys arise not only 
because of low-coordination and geometry differences but mainly because the symmetry 
enforces long Cu — Cu and Cu — Ag bonds, unlike the situation in Ag27Cuj. 

In Ag27Cuj, the charge density along Ag — Cu bonds whose length is ~ 2.6 A is even 
larger than that around Cil — Cu bonds, and certainly larger than that bridging Ag — Ag 
atoms. Nevertheless, Ag — Cu bonds whose length is of the order of that in bulk Ag^Cu, 
or even CusAg, are surrounded by an appreciably quenched charge density, explaining why 
the DOS of Cu atoms show low coordination features. 

We infer a hierarchy of bond strength intrinsic to Cu and Ag: {Cu — Cu > Cu — Ag > 
Ag — Ag}. The nanoalloy structures that enforce such a hierarchy also obey a bond-length 
order, {Ag — Ag > Cu — Ag > Cu — Cu}. Clearly, the strength and length of the bonds 
are correlated properties in any particular structure. Hence, for a given composition of a 
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nanoalloy, the global-minimum structure will be that which simultaneously maximizes the 
number of the strongest bonds while maintaining the bond-length order, and thus fulfilling 
the bond strength hierarchy that ultimately minimizes the energy. 
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TABLE I: This table contains six sets of 2-column sub-tables showing the distance from each type 
of atom in the nanoalloy to all its neighbors. In the right column of each sub-table appears the 
type of neighbor which is being referred and the number of such equivalent atoms at the same 
distance is shown in parenthesis. Notice that here equivalent atoms are not considered being so if 
they do not belong to the same layer. 



CuO AgO Cul Ag2 AgS AgA 



Type 




Type 




Type 




Type 




Type 




Type 




(NX) 


f/(A) 


(NN) 


f/(A) 


(NN) 


r/(A) 


(NN) 


r/(A) 


(NN) 


rf(A) 


(XX) 


r/(A) 


Cul(^) 


2.548 


CuO^^) 


2.715 


CuO(^) 


2.548 


CuO(^^ 


2.587 


Cul(^) 


2.726 


Cul(^) 


2.590 


Cu-l'-^^ 


2.548 


Ag3^^^ 


2.882 


Cu-l^^^ 


2.576 


A 

Ag-2^^> 


2.867 


CuO^^^ 


2.786 


A 


2.932 


Cu0(2) 


2.584 


A ^^^^ 

Ag-3(^) 
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Ag3(2'' 
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AgO*-^^ 


2.882 
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4.457 
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3.002 


Ag3(^^ 
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A ^fO\ 
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3.002 


Ag2^2) 


2.902 
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4.991 
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5.166 
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4.360 


CuO(^^ 


4.462 


Ag3(2) 


2.948 


AgO^^) 


5.691 
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A 
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4.425 


Ag-3^^^ 


6.721 
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6.866 
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CuO^^) 


5.277 
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FIG. [T] (Color online) Three-dimensional picture of the six types of non-equivalent atoms. 
The vertical arrow represents the 5-fold rotation axis, (a) Five CuO atoms sitting on the 
mirror plane of the cluster, zq = 3.88 A. Note that the Zj-coordinates are given with respect 
to the reference frame used in the figure and will define later interlayer distances; (b) The 
five AgO atoms also sit on the mirror plane, zq = 3.88 A; (c) Two Cul atoms sitting on the 
5-fold rotation axis symmetrically located above and below the mirror plane, zi= 5.17 A 
and z_i=2.59 A. Notice that CuO-pentagon fits in the A^fO-pentagon; (d) Ten Ag2 atoms 
form two pentagons symmetrically located above and below the mirror plane, Z2 = 5.31 A 
and z_2 = 2.44 A; (e) Ten Ag3 form two pentagons symmetrically located above and below 
the mirror plane, Z3 = 6.24 A and z_3 = 1.52 A; (f) The two AgA atoms sit on the 5-fold 
rotation axis symmetrically located above and below the mirror plane, Z4 = 7.76 A and z_4 

= oA. 

FIG. |2] Top view of Ag27Cui nanoalloy, perpendicular to the mirror plane. The side 
length of each pentagon as aM„ (see text). 

FIG. [3] (a) Bond coordination for Cu atoms of Ag2-rCuj compared with Cu bulk; (b) 
bond coordination for Ag atoms of Ag27Cuy compared with Ag bulk. 

FIG. [4] Calculated phonon spectra of the hypothetical LI2 phase of (a) Ag^Cu and (b) 
CusAg bulk alloys. 

FIG. [5] Band structure of Ag27Cui showing that the bands do not depend on k (see 
Ref. and that the rf-bands are concentrated between -5.0 and -1.5 eV below the Fermi 
level (see text). 

FIG. [6] (Color online) (a) Total and projected electronic DOS of Ag27CuT. The later 
corresponds to the contribution of s and d atomic states between -7 and 4.2 eV from the 
Fermi level, Ep, which is set equal to 0. The s-contribution is negligible up to Ep. The 
HOMO-LUMO gap in the ground state, A = 0.77 eV, is highlighted in red; (b) comparison 
between the total DOS of Ag2jCuj and that of Ag and Cu bulk; (c) contribution from 
the core {Cu) and shell {Ag) atoms to the projected DOS of Ag27Cui; (d) DOS of CusAg 
showing the ^-contribution from each species to the total DOS of the alloy in the LI2 phase; 
(e) DOS of Ag^Cu showing the ci-contribution from each species the total DOS of the alloy 
in the LI2 phase. 

FIG. [7] (Color online) Electronic DOS (rf-band) of Cu and Ag atoms situated in various 
environments and bond length (/); (a) Ag atoms in Ag^Cu {I = 2.87 A); (b) Ag atoms in 
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Cu^Ag [l = 2.70 A); (c) Ag atoms in Ag2'jCu'j; (d) Ag atoms in A free standing Ag{lll) 
monolayer [l = 2.94 A); (e) Ag atoms in bulk {I = 2.94 A); (f) Cu atoms in Ag^Cu {I = 
2.87 A); (g) Cu atoms in Cu^Ag {I = 2.70 A); (h) Cu atoms in Ag27Cu^; (i) Cu atoms in 
A free standing Cu{Hl) monolayer (/ = 2.59 A); (j) Cu atoms in bulk (/ = 2.59 A). 

FIG. [s] (Color online) Electronic DOS (rf-band) of Cu and Ag atoms in pure bulk and LI2 
bulk Ag — Cu alloys with expanded/contracted /: (a) Ag atoms in 2.4 % expanded Ag^Cu 
with I = lAg^M — 2-94 A; (b) Ag atoms in 9.8 % contracted Ag^Cu with I = Icu^^h: = 2.59 A; 
(c) Ag atoms in 8.9 % expanded Cu^Ag with / = lAgtuik'i W ^9 atoms in 5.5 % contracted 
CusAg with / = d{CuO-Cul)= 2.55 A; (e) Ag atoms 8.0 % expanded bulk Ag; (f) Ag atoms 
in 11.9 % contracted bulk Ag with / = Icutuik'' (s) 2.4 % expanded Ag^Cu; (h) Cu 

in 9.8 % contracted Ag^Cu; (i) Cu in 8.9 % expanded Cu^Ag; (j) Cu in 5.5 % contracted 
CusAg; (k) Cu in 8.0 % expanded bulk Cu; (1) Cu in 1.6 % contracted bulk Cu. 

FIG. [9] (Color online) The PDOS of the six non-equivalent types of atoms in the Ag27Cu7 
nanoalloy is compared with the DOS of pure bulk Cu (dashed lines) and Ag (dotted lines). 
Ep is shifted to for all of them; (a) AgO; (b) Ag2; (c) Ag3; (d) AgA (e) CuO; (f) Cul. 



FIG. 10 2D Charge density plot at a plane that contains the 5-fold rotation axis of 
Ag27Cu-j nanoalloy and is, therefore, perpendicular to its mirror plane. Atoms labelled with 
black color are precisely centered on that plane. 



FIG. [TT] 2D Charge density plot at the mirror plane of Ag27Cui. Atoms labelled with 
black color are precisely centered on that plane. 



FIG. 12 2D Charge density plot at a plane that contains the positions where CuO and Ag3 
atoms (labeled with black color) are centered. Atoms labelled with grey color are depicted 
in the section but are not centered precisely on that plane. 

FIG. [13] 2D Charge density plot at a plane that contains the positions where Ag2 and Ag3 
atoms (labeled with black color) are centered. Atoms labelled with grey color are depicted 
in the section but are not centered precisely on that plane. 



FIG. 14 2D Charge density plot at a plane which is parallel to the mirror plane of Ag27Cu7. 
The pentagonal layer composed of Ag3 atoms (labeled with black color) is contained in this 
plane. Atoms labelled with grey color are depicted in the section but are not centered 
precisely on that plane. 



FIG. 15 2D Charge density plot at a plane that contains the positions where AgO and Ag2 
atoms (labeled with black color) are centered. Atoms labelled with grey color are depicted 



25 



in the section but are not centered precisely on that plane. 
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